Pseudo-gap pairing in ultracold Fermi atoms 
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The BEC-BCS crossover in ultracold Fermi gases créâtes an idéal environment to enrich our knowl- 
edge of many-body Systems. It is relevant to a wide range of fields from condensed matter to 
astrophysics. The nature of pairing in strongly interacting Fermi gases can be readily studied. This 
aids our understanding of related problems in high-T^ superconductors, whose mechanism is still 
under debate due to the large interaction parameter. Here, we calculate the dynamical properties of 
a normal, trapped strongly correlated Fermi gas, by developing a quantum cluster expansion. Our 
calculations for the single-particle spectral function agrée with récent rf spectroscopy measurements, 
and clearly demonstrate pseudogap pairing in the strongly interacting régime. 
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C/2 

I ' 



o 
o 



(N 
> 
00 

m 



m 
O 

o 



X 



Ultracold Fermi atom experimeiits allow inter-particle 
interactions, geometries, and spin species of atomic gases 
to be precisely controlled and tuned at will, leading to 
an exactly known model Hamiltonian ^j] . Since rapid ex- 
périmental progress is able to provide accurate data on 
both static [24j| and dynamic fs'-lq properties, this type 
of Quantum Simulation can definitively settle fundamen- 
tal issues. An interesting example is momentum-resolved 
rf spectroscopy at the BEC-BCS crossover @, which 
can reveal single-particle excitation gaps or pseudogaps 
in the normal state, as in analogous high- Te supercon- 
ductor Systems |9l4l3j. 

By contrast, theoretical progress [Ï^JÏs'l in describing 
strong pairing fluctuations at the BEC-BCS crossover is 
notoriously difficult due to the lack of a small interac- 
tion parameter fîÔ'l. The situation is most severe for 
dynamical properties, where différent crossover théories 
lead to qualitatively différent prédictions. For the single- 
particle spectral function in the strong-coupling régime, 
some crossover théories predict a pseudogap - the pre- 
cursor of fermionic pairing in the normal state above Te 
23"24| - while some others (251 claim no such efFects. 
Quantum Monte Carlo simulation of the spectral func- 
tion is not conclusive |26| . 

In this Letter, we develop a quantum virial or clus- 
ter expansion to solve this délicate problem of dynamic 
properties for a normal, trapped, and strongly interact- 
ing Fermi gas. The advantages of this method are clear. 
First, the expansion has a controUable parameter [27l . [28| : 
the fugacity z = exp(^/fcsT) is small at high tempéra- 
tures T. Second, the expansion coefficient or function at 
the n-th order (n ^ 2) is entirely determined by knowl- 
edge of a n-particle cluster. Multiparticle corrélations, 
which are missing in most current crossover théories, 
can be easily accounted for and improved. Finally, the 
method is easily capable of treating external potentials 
or traps. We have estimated the validity of the expan- 
sion by comparing its prédictions with expérimental data 



for thermodynamics [19|]. At the cusp of the crossover, 
where the s-wave scattering length Og diverges (unitarity 
liniit [2^), it is applicable down to 0.4Tf for a trapped 
Fermi gas [19.] . 

To develop the cluster expansion for gênerai dynamical 
properties, let us consider two arbitrary linear operators 
of physical interest, R and S, and the related Green func- 
tion or corrélation function at différent space-time points. 



G(r,r';r) = - 



Tr 



^-m-f^J^)R{j.^r) S+ (r') 



(1) 



where at finite températures we are working with an 
imaginary time r in the interval < r < /3 = l/fc^T. 
At high températures, both numerator and denomina- 
tor may be expanded into the powers of z ^ 1, lead- 
ing to G{v,v';t) = (Xo + zXi + + zQi + 
•••) = Xq + ziXi- XqQi) + ■■■ , where X„ = - 
Tr„[e~^'^iî (r,T) ^^(r')] is the expansion function and 
Qn =Tr„[e~''^] is the cluster partition function. The 
above expansion is to be referred to as the cluster expan- 
sion of corrélation function, G (r, r'; t) = G'-*^' (r, r'; r) -|- 
zG'^^^ (r, r'; r) -I- • • • , where. 



G(°)(r,r';r) = Xq, 
GW(r,r';r) = Xi-XqQi, 



etc. 



(2) 



The experimentally measured spectral function v4(k,cj) 
can be calculated from the corrélation function via an- 
alytic continuation, so that we may write accordingly, 
A (k, uj) = (k, u) + (k, cj) -}-••• . The calcu- 
lation of the n-th expansion coefficient G*-"-' (r,r';T) or 
(k,w) requires the knowledge of solutions up to the 
ri-body problem, including both energy levels and wave- 
functions (see Appendix). In this work we aim to calcu- 
late the leading efFects, which give the 2"'^-order expan- 
sion function. The next-order expansion function [281], is 
not treated here for simplicity. 

To calculate the single-particle spectral function 
A{k, Lj), we take the annihilation field operators '^a{Y, t) 
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Figure 1: (Color online) Contour plots of the occupied 
spectral intensity at crossover. The intensity I{uj) = 
A(k,uj)f{uj)k'^/{2iT^) increases from blue {W~^Imax) to red 
(Imax) in a logarithmic scale. The calculations were per- 
formed with harmonie traps at T = 0.7Tp and l/{kFas) = 
+ 1, 0, —1, with a resulting fugacity at the trap center of 
z ~ 0.14, 0.42, and 0.48, respectively. 

{a =t,i) for R 01 S (see Appendix). Fig. 1 shows con- 
tour plots of the occupied spectral intensity of a trapped 
Fermi gas in the crossover at T = Q.TTp. At this tem- 
pérature, our results are quantitatively reliable [ï^ . We 
observe that, in addition to the response from cohérent 
Landau quasiparticles (black lines) , there is a broad inco- 
hérent spectral weight centered about lo + ^ ~ — Ck — es 
(white dashed lines), where ek — fi.^fc^/(2m) and es = 

/{maj.) is the binding energy. Thus, the spectra clearly 
exhibit a gap-like double peak structure in the normal 
state. This is a remarkable feature: the dispersion at 
négative énergies seems to follow the BCS-like dispersion 
curve, uj = —a/ (fk — m)^ + A^, and behaves as if the gas 
was superconducting, even though we are above the crit- 
ical température Te. Therefore, the incohérent spectral 
weight indicates the tendency of pseudogap: the precur- 
sor of fermionic pairing due to strong attractions, i.e., it 
arises from the atoms in the paired state or "molécules". 
The pairing response is very broad in energy and bends 
down towards lower energy for increasing k. At large 
\|{kF<ls)^ the width is of order -^/maxIfesT, êfIêIc- 

The incohérent spectral weight found by our leading 
cluster expansion is a universal feature of interacting 
Fermi gases. At large momentum fc 3> /cf; it is related to 
the universal 1/fc^ tail of momentum distribution [30, 31], 
rio.(k) = j dio A{k, uj) f I/fc^, arising from the 
short-range two-body corrélations. A nonzero contact T 
therefore necessarily implies a finite spectral weight at 
négative énergies. At fc ^ fci?, we anticipate that many- 
body corrélations will become increasingly important. 

For an absolute comparison with experiment Q, we 
perform calculations using realistic expérimental param- 
eters, including the measurement resolution. Fig. 2 
présents the results on the BEC side of crossover with 
l/ikpas) = 1.1. The température T = 0.45rF is esti- 
mated from an initial température Ti = O.lTTp obtained 
before the field sweep to the BEC side The exper- 
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Figure 2: (Color online) Single-particle excitation spectra on 
the BEC side of crossover. a-c, Cluster expansion prédictions 
(2 ~ 0.1 and ^ ~ — LOSef). d-e, Corresponding expérimen- 
tal data 3. a, The linear-scale intensity map. Our results 
were convoluted with a Gaussian broadening curve of width 
a — 0.22eF, to account for the measurement resolution 
The black line shows upper free-atom dispersion. The red 
dashed line is the lower dispersion curve of molécules, ob- 
tained via fitting each fixed-fc energy distribution curve (in 
b) with a two Gaussian distribution. It agrée fairly well with 
the expérimental resuit (white symbols). b, Energy distribu- 
tion curves for selected values of k. c, The occupied density 
of state (DOS). Blue dashed lines show the expérimental peak 
positions. 



imentally observed upper and lower features, caused re- 
spectively by unpaired atoms and molécules, are faith- 
fully reproduced. In particular, the expérimental data 
for the quasiparticle dispersion of molécules, marked by 
white symbols, agrées with our theory (lower red dashed 
line) . There is also a qualitative agreement for the energy 
distribution curves (Figs. 2b and 2e) and the occupied 
density of states (Figs. 2c and 2f). A narrow peak due to 
free atoms and a broader feature due to molécules are re- 
produced theoretically with very similar width at nearly 
the same position. It is impressive that our simple quan- 
tum cluster expansion is able to capture the main feature 
of the expérimental spectra. In contrast, a more compli- 
cated crossover theory with adjustable parameters fails 
to account for the free-atom contribution with similar 
parameters [24|. 
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Figure 3: (Color online) Single-particle excitation spectra of 
a strongly interacting Fermi gas. a-c, Cluster expansion pré- 
dictions (z ~ 6 and fi ~ 0.37eF). d-e, Corresponding expéri- 
mental data Q. In e, for the expérimental energy distribution 
curves, we use a larger value of k (i.e., enlarged by a factor of 
5/3) to account for a scaling discrepancy due to many-body 
corrélations. 



Fig. 3 reports the spectra in the unitarity hmit at 
the critical température Te ~ 0.2Tf. At such low tem- 
pératures, the use of a cluster expansion becomes highly 
questionable as the fugacity at the center z ~ 6 ^ 1. 
Nevertheless, we find that the dispersion curve is low- 
ered by the attractions by an amount comparable to the 
Fermi energy , as shown clearly by the red dashed line 
in Fig. 3a. The calculated energy distribution curves 
bifurcates from a single peak with increasing k and be- 
comes dominated by the lower molecular branch (Fig. 
3b), which eventually leads to the bending back of the 
dispersion curve to négative energy. This picture is sug- 
gestive of the existence of a pseudogap and is consistent 
with the expérimental findings (Fig. 3e). This surpris- 
ingly good agreement merits further investigation. We 
conjecture that even at thèse relatively low températures 
the virial expansion captures the dominant two-body cor- 
relations measured in thèse experiments, apart from a 
possible overall scaling factor due to the missing higher- 
order terms. 

To study the température dependence of pseudogap 
pairing, we show in Fig. 4 the unitary spectral inten- 
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Figure 4: (Color online) Température dependence of the spec- 
tral intensity at unitarity. a-c, The intensity increases from 
blue [10^'^ Imax) to red (Imax) in a logarithmic scale. d and 
e. Température dependence of the fugacity and chemical po- 
tential at unitarity. Our leading cluster expansion appears to 
be applicable at T > OATf, where 2; ~ 1 or ^ ~ |19l |. 



sity at T/Tc — 1.2, 1.5 and 2.0, to be compared with 
the most récent measurement at JILA [8|. In this loga- 
rithmic scale, a pseudogap structure is clearly visible at 
about k ^ 0.5kp and w + /i ^ —O.bep. The response is 
reduced with increasing température and becomes signif- 
icantly weak once T > 2Tc. However, it seems difficult to 
accurately détermine a characteristic température above 
which the pseudogap response disappears. 

In conclusion, the results obtained here provide a good 
qualitative explanation for the récent single-particle spec- 
tral function measurement in trapped strongly interact- 
ing ultracold fermions. It demonstrates the precursor of 
fermionic pairing in the normal but strongly interacting 
régime. In the near future, tomographic rf spectroscopy 
may be used to reveal locally the homogeneous spectral 
function. Our calculation can be extended to a uniform 
gas, for which we expect a stronger pairing response and a 
wider température window for pseudogap. Our method is 
also applicable directly to multi-component atomic Fermi 
gases and can be used to understand the intriguing triplet 
and quadruplet pairing in their dynamic responses. 
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Appendix. General rules of cluster expansion. — To 
calculate the Green function or corrélation function, it is 
convenient to separate out the contribution arising from 
interactions. To this aim, for any physical quantity Q 
we may write Q = {S}^^^ + Q^^\ where the superscript 
"TV" in Q(^) dénotes the part of a non-interacting System 
having the same fugacity. The operator {\^^^ thcn picks 
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up the residues due to interactions. We then may write, 
G (r, r'; r) = {G (r, r'; t)^^ + G^^) (r, r'; r) , 

where {G}^^^ can be expanded in terms of 

For the single-particle spectral function, we détermine 
first the Green's function Gg-cr' (r, r'; t) by taking R = ^„ 
and S = ^ct' • Then, we take the Fourier transformation 
with respect to r, to obtain Gaa' (r, r'; i'^n)- The spectral 
function is calculated via analytic continuation, 

Acra' (r,r'; w) = --ImGaa' (r,r' uj + iQ+) . 

A final Fourier transform on r — r' leads to A^^i (k, , as 
measured experimentally. For a normal, balanced Fermi 
gas, A^^ ^ = A{k, w) and A^^ = 0. 

Leading expansion of spectral function. 
term of {G^^{r,r';T)}'^^^ takes the form. 



Tn 



-13-HtH 



*t (r) e-^'^^ (r')] } 



{I) 



The trace has to be taken over ail the single-particle 
States (i.e., with energy e^) for a spin-down fermion. 
We insert in the bracket an identity \Q) {Q\ = 1, 
where Q refers to the "paired" state (i.e., $q with energy 
Eq) for two fermions with unlike spins. It is straightfor- 
ward to show that, at the leading order, 



il) 



-ze' 



-f!cp+T{ep~EQ) p 



pQ 



where FpQ = J dridr2ip*{ri)^Q (r,ri)$^ (r', rj) i/'p(r2). 
Accordingly, the leading interaction correction to the 
spectral function, {A(k, is given by. 



(1 + ^ 



E 



s {io + Ép - Eq + fi) e 



F ru 



il) 



where FpQ(k) = J drdrie-'^ ''i^*{ri)<èQ (r, ri). 

In an isotropic harmonie trap with frequency ujq, we 
solve exactly the two-fermion problem for relative wave- 
functions [28| and obtain {A (k, uj)}^-^'> using the above 
procédure. In the end, we calculate 



/(k,c.) 



27r2 



as measured experimentally [6, ^8]. Hère, / (i^) — 
l/{e^'^ + 1) and the idéal spectral function A'-^'> — 
4-\/27r/ (m'^/^Wg) {uj + fi — ek)^^^- To account for the ex- 
périmental resolution, we further convolute /(k, w) with 
a Gaussian broadening curve. 

In the BEC limit, we may show analytically that. 



{A}'^^^ f (X exp ~I3 (V^k - w - /.j - e_B - v^) 



where eb = 1^? /{mal) ^he binding energy. Thus, at 
large k the intensity peaks at w -|- = — ek — es, with a 
width ^ ^/ksTïû- 
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